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Abstract. We present the first high resolution X-ray 
spectrum of the bright 04Ief supergiant star C, Puppis, 
obtained with the Reflection Grating Spectrometer on- 
board XMM-Newton . The spectrum exhibits bright emis- 
sion lines of hydrogen-like and helium-like ions of nitrogen, 
oxygen, neon, magnesium, and silicon, as well as neon-like 
ions of iron. The lines are all significantly resolved, with 
characteristic velocity widths of order 1000 — 1500 km s 1 . 
The nitrogen lines are especially strong, and indicate that 
the shocked gas in the wind is mixed with CNO-burned 
material, as has been previously inferred for the atmo- 
sphere of this star from ultraviolet spectra. We find that 
the forbidden to intercombination line ratios within the 
helium-like triplets are anomalously low for N VI, 0 VII, 
and Ne IX. While this is sometimes indicative of high elec- 
tron density, we show that in this case, it is instead caused 
by the intense ultraviolet radiation field of the star. We 
use this interpretation to derive constraints on the loca- 
tion of the X-ray emitting shocks within the wind that 
agree remarkably well with current theoretical models for 
this system. 
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1. Introduction 

The initial discovery of X-ray emission from O stars with 
the Einstein Observatory in the late 1970s (Harnden et al. 
1979), sparked a vigorous field of research aimed at bet- 
ter understanding the production of hot gas in such sys- 
tems. Most current models invoke hydrodynamic shocks 
resulting from intrinsic instabilities in the massive, radia- 
tively driven winds from these stars (for a detailed review 
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see Feldmeier et al. 1997). Support for this picture comes 
from the fact that the X-ray flux is not highly absorbed 
at low energies (Cassinelli & Swank 1983; Corcoran et al. 
1993), which is expected if the X-ray emitting gas is dis- 
tributed throughout the wind rather than in some form of 
hot corona in the outer atmosphere of the star. However, 
the X-ray observations to date have not been especially 
constraining for stellar wind models. This is primarily due 
to the low spectral resolution of the available nondisper- 
sive detectors, which has precluded the study of individual 
atomic features so crucial to the unambiguous determina- 
tion of physical conditions in the shocked gas. 

In this Letter, we present one of the first high reso- 
lution X-ray spectra of an early-type star, the 04Ief su- 
pergiant £ Puppis, which we have obtained with the Re- 
flection Grating Spectrometer (RGS) experiment on the 
XMM-Newton Observatory. ( Pup is an excellent target 
for such work since it is the brightest O star in the sky, and 
has consequently been very well studied at longer wave- 
lengths. In particular, Pauldrach et al. (1994) have devel- 
oped a very detailed NLTE model of the atmosphere and 
wind of this star, constrained by high quality ultraviolet 
spectra from Copernicus and RJE. They find a mass loss 
rate of ~ 5.1 x 10“ 6 M e yr -1 using an effective tempera- 
ture ~ 42 000 K. The derived abundances indicate that the 
atmosphere is mixed with CNO-processed material, con- 
sistent with the theoretical picture of a highly evolved star 
near the end of core hydrogen burning. Hillier et al. (1993) 
found that the Pauldrach et al. model is consistent with 
the ROSAT X-ray spectrum of C Pup, but only if the X- 
ray emission arises in shocks distributed throughout the 
wind. At the very lowest energies (< 200 eV), they find 
that the emitting material must be > 100R* away from 
the star. This is a consequence of the fact that in the Paul- 
drach et al. model, helium is mostly only singly ionized in 
the outer regions of the wind, so that the photoelectric 
opacity is very high at soft X-ray energies. 
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Our XMM-Newton RGS spectrum is dominated by 
broad emission lines of mostly hydrogen-like and helium- 
like charge states of nitrogen, oxygen, neon, magnesium, 
and silicon, and neon-like ions of iron. The data provide 
a number of important constraints on the nature and lo- 
cation of the X-ray emitting material in the ( Pup wind. 
Of particular interest is the fact that we see a suppression 
of the forbidden line and enhancement of the intercombi- 
nation line in the helium-like triplets of nitrogen, oxygen, 
and neon. We show that this is a natural consequence 
of the intense ultraviolet radiation field in the wind, and 
that it allows us to place constraints on the location of 
the X-ray emitting shocks relative to the star. Additional 
constraints come from the emission line velocity profiles. 
We show that the data agree well with the predictions of 
the Pauldrach et al. model. 

In section 2, we describe the details of the RGS obser- 
vation, the nature of our data reduction and analysis, and 
the key observational features of the spectrum. In section 
3, we consider the implications of our results for our un- 
derstanding of C Pup, and of the X-ray emission from O 
stars in general. 


2. Observations and Data Analysis 

The RGS covers the wavelength range of 5 to 35 A with 
a resolution of 0.05 A, and a peak effective area of about 
140cm 2 at 15 A (see den Herder et al. 2000; this vol.). 
C Pup was observed for 57.4 ks on 2000 June 8. The data 
were processed with the XMM-Newton Science Analysis 
Software (SAS). Filters were applied in dispersion channel 
versus CCD pulse height space to separate the spectral or- 
ders, and the source region was separated with a 1' spatial 
filter. The background spectrum was obtained by taking 
events from a region spatially offset from the source. The 
wavelengths assigned to the dispersion channels are based 
on the pointing and geometry of the telescope and are ac- 
curate to ~ 0.008 A. The effective area was simulated with 
a Monte Carlo technique using the response matrix of the 
instrument and the exposure maps produced by the SAS. 
Based on ground calibration, we expect the uncertainty 
in the effective area to be less than 10% above 9 A and 
at most 20% for shorter wavelengths. A fluxed spectrum, 
corrected for effective area, with the two first order spec- 
tra added together to maximize statistics, is presented in 
Fig. 1. The spectrum was also extracted for analysis in 
XSPEC using standard SAS routines. 

Due to a large solar flare event that occurred close to 
the time of the £ Pup observation, all three EPIC detectors 
and the Optical Monitor on XMM-Newton were switched 
off while this source was observed. Thus, only data from 
RGS are available for analysis. 
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Fig. 1. First order background subtracted spectrum. It has 
been corrected for the effective area of the instrument. 
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Fig. 2. Light curve of first order events in the range 6.2 to 13.8 
A, corresponding to 0.9 to 2.0 keV. Events from both instru- 
ments have been binned together. The time dependence is plot- 
ted as a function of the 16.67 hr period reported in Berghofer 
et al. 



2.1. Light Curve 

Previous ROSAT obervations of £ Pup (Berghofer et 
al. 1996) show evidence for a 16.67 h period of variability 
in the X-ray emission between 0.9 and 2.0 keV. We have 
extracted a light curve (Fig. 2) of the events in this energy 
range in the RGS data set. There is no significant observed 
variation on this time scale apparent to the eye, and the 
fit to a constant intensity yields xt — 0.9. However, our 
upper limit to the percentage variation at this period is 
not inconsistent with the Berghofer et al. detection. 

2.2. Emission Line Intensities and Emission Measure 
Analysis 

The spectrum displayed in Fig. 1 is composed almost en- 
tirely of emission lines, with a weak underlying contin- 
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uum. The lines are all broadened, with characteristic ve- 
locity widths of order 1000 km s' 1 . In nitrogen, oxygen, 
and neon the He-like forbidden lines are very weak and 
the intercombination lines are bright. The line fluxes (Ta- 
ble 1) were measured by taking the integrated flux of the 
spectrum over the line and subtracting a corresponding 
amount of continuum flux. Lines that originate from the 
same ion and that are blended were evaluated as one com- 
plex (for example, the He-like triplets or the Fe XVII emis- 
sion around 15 or 17 A). The continuum strength was de- 
termined by taking the flux of a spectral region free of lines 
but near the line in question. In cases where this was not 
possible, the continuum strength was interpolated from 
the strength in other regions of the spectrum. Some of the 
measurements were complicated by the presence of over- 
lapping lines. When possible, the flux of these other lines 
was estimated from the flux of lines originating from the 
same ion by comparing with line power ratios. 

We calculated the emission measure (see Fig. 3) for 
each ion assuming solar abundances (Anders & Grevesse 
1989) and a temperature given by the temperature of for- 
mation for the dominant lines from that charge state. The 
emission measure is 


EM = F Unc 


4 ncP 

F]ineA/i 


( 1 ) 


where F| in e is the observed flux in the line, d is the dis- 
tance to the source, Pane is the line power, A is the ele- 
mental abundance, and f\ is the ion fraction evaluated at 
the temperature of formation. We used line powers from 
the APEC code (Smith & Brickhouse 2000), which in- 
cludes ion fractions from Mazzotta et al. (1998). We take 
d = 450 pc (Schaerer et al. 1997). The line fluxes must 
be corrected for interstellar absorption, although this is 
a minor effect. We take Nh = 10 2O cm" 2 (Chlebowski et 
al. 1989) and we use cross sections from Morrison & Mc- 



Fig. 3. Emission measure calculated from each line as a func- 
tion of Tform- The horizontal error bars indicate the tempera- 
ture range over which Pii„ e /» is at least half its maximum value. 
Solar abundances are assumed. 


Table 1. Measured fluxes for prominent emission lines with 
rest wavelengths in A. 


Line 

Flux 

(photons cm - 2 s~ 1 ) 

Si XIII (6.65,6.69,6.74) 

(1.1 ±0.4) x 10 _,: 

Mg XII (8.42) 

(4.1 ±2.8) x 10 -3 

Mg XI (9.17,9.23,9.31) 

(2.0 ±0.7) x 10 -2 

Ne X (12.13) 

(1.8 ±0.7) x 10" 2 

Ne IX (13.45, 13.55, 13.70) 

(5.0 ± 1.0) x 10“ 2 

Fe XVII (15.01, 15.26) 

(7.4 ±0.7) x 10“ 2 

Fe XVII (16.78, 17.05, 17.10) 

(4.7 ± 0.8) x 10~ 2 

O VIII (18.97) 

(3.50 ±0.22) x 10 -2 

O VII (21.60,21.80,22.10) 

(5.1 ± 0.4) x 10 -2 

N VII (24.78) 

(6.4 ± 0.5) x 10 -2 

N VI (28.78,29.08,29.53) 

(1.13 ±0.11) x 10 _1 

C VI (33.74) 

(6.4 ± 2.8) x 10 -3 


Cammon (1983). 

We find that the emission measures derived from the 
nitrogen emission lines are at least an order of magnitude 
greater than those for carbon and oxygen, which both 
have temperature of formation ranges that overlap with 
that of nitrogen. This indicates that the ratios A^/Ao 
and An /Ac are substantially higher than solar, even al- 
lowing for a factor of two uncertainty due to the crudeness 
of the emission measure analysis. This result is consistent 
with the inference of atmospheric abundances by Paul- 
drach et al. (1994), based on their analysis of UV spectra. 
It indicates that the wind material in £ Pup has been 
significantly mixed with matter that has undergone CNO 
burning in the stellar interior. 

The emission measure expected for a smooth, spher- 
ically symmetric wind with parameters appropriate to 
C Pup is EM = 6.5 x 10 6O cm" 3 . The fact that we find 
EMs that are lower by 2 to 3 orders of magnitude implies 
that only a small fraction of the wind material is heated 
to X-ray emitting temperatures. 


2.3. Continuum emission analysis 

The continuum emission has a very low intensity relative 
to the bright line emission. To ensure that the contin- 
uum that we see in the spectrum is real, and not, for 
example, an artifact of faulty background subtraction, we 
plotted the first order events in a region free of emission 
lines as a function of their spatial distribution in the cross- 
dispersion direction, without subtracting the background. 
The peak at the location of the source spectrum is clearly 
visible, and indicates the presence of true continuum, well 
above background. The detected continuum also signifi- 
cantly exceeds the scattered light contribution from the 
ensemble of detected emission lines. 

We looked for evidence of discrete photoelectric ab- 
sorption edges in the continuum. Although there axe no 
edges evident, the constraints are weak due to the low 
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Table 2. Upper limits on the strengths of the K-edges of Ne, 
O, and N. 


Edge 

Wavelength 

(A) 

Upper limit on to 

Ne X 

9.10 

0.5 

Ne IX 

10.37 

0.5 

O VIII 

14.23 

0.2 

O VII 

16.78 

0.3 

N VII 

18.59 

0.5 

N VI 

22.46 

0.5 


strength of the continuum emission. The upper limits on 
the optical depths are given in Table 2. These upper lim- 
its are incompatible with the detection of a strong edge 
feature near 0.6 keV reported by Corcoran et al. (1993), 
but that is perhaps not surprising, considering the com- 
plexity of the spectrum and the low spectral resolution of 
their measurements. 

The continuum is so weak that it is difficult to ascer- 
tain its overall shape, especially in the region from 8 - 
17 A where there axe many lines. We tried fitting with 
a bremsstrahlung model in XSPEC, but the results were 
inconclusive. 


2.4- He-like triplet ratios 

The forbidden to intercombination line ratios, R = //t, 
for helium-like oxygen and neon were obtained by fitting 
in XSPEC. In each case, we fit a gaussian to the Lya 
line of the respective element, and used that profile for 
each of the three components in the He-like triplet. We 
also take into account both the low level continuum and 
lower flux lines in the near vicinity. For Ne IX we find R = 
0.34 ±0.11 and for O VII we find R = 0.19 ±0.08. The fit 
to Ne IX is shown in Fig. 4. These values are well below 
the expected values for low density plasmas in collisional 
equilibrium (see Mewe et al. 2000). A similar analysis did 
not work very well for the N VI He-like triplet because of 
the complexity of both the Lya and He-like profiles for 
that element. Nevertheless, it is clear from the data that 
the forbidden line is strongly suppressed for N VI as well. 
The Si XIII and Mg XI triplets are too blended to allow 
the forbidden line and intercombination line to be quan- 
titatively separated. However, there is clearly significant 
emission in Mg XI from the forbidden transition, at a level 
that is not incompatible with the low density limit. 

The conversion of forbidden line to intercombination 
line emission in high density plasmas is a well known ef- 
fect. However, this can also occur at much lower densi- 
ties if the plasma is exposed to a strong UV radiation 
field (Mewe & Schrijver 1978). To produce a low R ra- 
tio, electrons populating the 2 3 *Si state must be excited 
into the 2 3 P state. The state 2 3 S X is metastable, but 





Fig. 4. The Ne IX triplet with best fit model. Only RGS 2 
data is shown for clarity, but both were used in the fit. There 
is also a Fe XVII line at 13.825 A. The rest wavelengths of the 
three lines are at 13.45 A (r), 13.55 A (i), and 13.70 A (f). 


it is easily populated at collisional equilibrium temper- 
atures, and the corresponding emission line intensity is 
normally strong. When the excitation rate to the 2 3 P lev- 
els becomes comparable to the decay rate to the ground 
state, the forbidden line is suppressed. In a high density 
plasma, this occurs because the collisional excitation rate 
from 2 3 Si to 2 3 P competes effectively with radiative de- 
cay. However, since this is a dipole transition, it can also 
be photoexcited if the ambient U V flux at the appropriate 
wavelength is sufficiently high. To calculate the photoex- 
citation rate, we estimated the emission from £ Pup for 
the frequencies of the 2 3 S X to 2 3 P transitions for N, 
O, Ne, Mg, and Si. We assumed a blackbody spectrum 
with ~ 42000K (Pauldrach et al. 1994). The rate 
of photoionization is given by Ppe = PV~~r/> where / is 
the oscillator strength. We used oscillator strengths from 
Cann and Thakkar (1992) and Sanders and Knight (1989). 

The flux is given by F v = 2ir(l - yj\ - ) 2 )J„, where 

I v is the specific intensity for the blackbody. 

In Table 3 we list the decay rates of 2 Z S\ and the 
photoexcitation rates from 2 3 S x to 2 3 P evaluated at the 
photosphere. It is clear from this calculation that the for- 
bidden line suppression we observe is in fact due to pho- 
toexcitation from < Pup’s high UV flux, and not due to 
collisional excitation at high densities, as long as the emit- 
ting regions are close enough that the U V is not sufficiently 
diluted. We can thus place constraints on the location of 
shock formation for each of the observed lines. Since the 
forbidden lines in N, O, and Ne are strongly suppressed, 
the emission must occur at radii smaller than the radii at 
which the photoexcitation and decay rates are equal. The 
emission from Mg and Si must occur at larger radii than 
the values given for those elements in the table. 


I 
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Table 3. Comparison of photoexcitation and decay rates of 
the 2 3 5 state. Decay rates are from Drake (1971). 


Ion 

Decay Rate 
2 3 S — ► 1 *5 

(O 

Photoexcitation 
Rate (R = R*) 
2 3 S — ► 2 3 P 
(*-') 

Radius 
at which 
R PE = -Rdecay 

(R.) 

Si XIII 

3.56 x 10 b 

4.96 x 10 s 

1.04 

Mg XI 

7.24 x 10 4 

7.53 x 10 5 

2.3 

Ne IX 

1.09 x 10 4 

1.13 x 10 6 

7.2 

O VII 

1.04 x 10 3 

1.69 x 10® 

29 

N VI 

2.53 x 10 2 

2.03 x 10 6 

63 


2.5 . Line profile analysis 

The projected velocity (v p ) profile for a thin spherical shell 
(with a single radial velocity) is flat, or = (const.). 
We expect the emission line profile to appear as a con- 
volution of the radial emission intensity with this flat 
projected velocity profile. The velocity of each shell is 
given as a function of radius by the conventional /3-model: 
v(r) «Uoo(l -r 0 /r) p (Lamers & Cassinelli 1999) where 
/3 « 0.8 and r 0 « R+. Lines formed at larger radii will 
therefore appear broader than lines formed close to the 
star. Furthermore, lines originating from larger radii are 
formed over a region with a small velocity gradient. Ther- 
fore, we expect these lines to appear more flat- topped than 
lines originating closer to the star. This is apparent in the 
observed Lya lines, which are plotted in velocity space in 
Fig. 5. The N VII peak is noticeably broader and has a 
substantially different shape than the other peaks in the 
plot. It also shows evidence for resolved, discrete struc- 
ture, given the resolution of the instrument. Note that 
other lines overlap with the Ne X line; it does not have a 
red shoulder. 

3. Discussion 

The RGS spectrum we have presented provides very 

strong confirmation for a number of key aspects of conven- 
tional models of X-ray emission from early-type stars in 
general, and for C Pup in particular. The emission is dom- 
inated by broad emission lines indicative of hot plasma 
that is flowing outward with the wind. We do not see ev- 
idence for strong attenuation at low energies, which con- 
firms that the X-ray emitting regions are not concentrated 
at the base of the wind, but instead are distributed out to 
large radii. These very simple results agree well with the 
predictions of wind instability models. 

As we have shown, the suppression of the forbidden 
lines in the He-like triplets of low Z elements allows us to 
derive upper limits to the radii of the emitting shocks in 
each case. It is interesting to compare these values with 
Fig. 1 of Hillier et al. (1993), which is a plot of the radius 
at which optical depth unity is reached in the wind as a 



Fig. 5. Ly a lines in velocity space. Nitrogen is the solid line, 
oxygen is dashed, and neon is dotted. The intensities have been 
renormalized for comparative purposes. The error bar is rep- 
resentative for the nitrogen line. Note the discrete structure 
in the nitrogen peak. The shoulder on the Ne X line is not a 
velocity feature; it is an emission line from Fe XVII. 


function of X-ray energy for the NLTE model of C Pup 
calculated by Pauldrach et al. (1994). At the energies of 
the N VI, O VII, Ne IX, Mg XI, and Si XIII lines, unit 
optical depth is achieved at 22, 23, 9, 3.5, and 2.5 stellar 
radii, respectively, in this model. These values are quite 
compatible with both our derived upper limits for N VI, O 
VII, and Ne IX, and our derived lower limits for Mg XI and 
Si XIII. Since the density in the wind drops off like r~ 2 , 
and the emissivity is proportional to n 2 , it is perhaps not 
surprising that we infer radii characteristic of the smallest 
radius at which the overlying wind is still transparent to 
the respective line, in each case. Thus our results are in 
excellent agreement with the Pauldrach et al. predictions. 

Further support for this picture comes from the ob- 
served velocity profiles. The higher Z lines have character- 
istic widths ~ 1000 km s“ l , whereas the N VII line is dis- 
tinctly broader. The terminal velocity in the £ Pup wind 
is 2260kms _1 (Groenewegen et al. 1989), so the higher Z 
lines are most likely emitted at only a few stellar radii, 
whereas N VII can come from considerably further out. It 
is also noteworthy that the N VII line is much less smooth 
than the higher Z lines. One might expect the shock struc- 
ture to be far more chaotic in the outer regions of the wind, 
leading to complex line profiles. 

We have shown that the respective line intensities 
are suggestive of significant enhancements of the nitrogen 
abundance relative to carbon and oxygen, as one would ex- 
pect for CNO-processed material. This, again, agrees well 
with the Pauldrach et al. model. Meynet & Maeder (2000) 
recently presented new evolutionary models for rotating 
single stars. They found that rotational mixing produces 
a significant surface helium and nitrogen enhancement. 
Meynet & Maeder suggested that stars with enhanced He- 
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abundances and large projected rotational velocities are 
natural descendants of very fast main-sequence rotating 
stars. In these stars, the chemical enrichment at the sur- 
face is very fast and as a result of the strong rotational 
mixing the chemical structure of these stars could be near 
homogeneity. With its large projected rotational velocity 
of v sin i ~ 203 km s" 1 (Penny 1996), (Pup most probably 
falls into this category. 
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